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Introduction 

A digital  signal  processing  laboratory  has  recently  begun  operation 
in  the  Electronics  Section,  Space  Technology  Branch,  Space  Systems 
Division.  This  report  details  the  general  design  approach  of  hardware  and 
software  systems  and  the  present  operating  capability.  Block  diagrams 
of  the  hardware  and  software  systems  are  provided  along  with  an  example 
of  a recent  task . 

Hardware  Systems  Design 

The  digital  signal  processing  laboratory  (Fig.  1,  2)  is  built  around 
a Varian  V-73  disk  based  mini-computer  system.1  This  computer  is  a 16-bit 
machine  with  dual  asynchronous  Input/Output  (I/O)  busses  and  dual  ported 
memory.  It  has  I/O  data  transfer  rates  of  up  to  2.9  million  words  per  second. 
Another  feature  of  this  machine  is  a 512  64-bit  word  Writable  Control  Store 
(WCS).  The  Disk  has  a 1.7  million  16-bit  word  capacity,  expandable  to  8 
million  words.  The  present  computer  memory  size  is  16K  words. 

Other  key  elements  of  the  laboratory  are:  a 5 MHz  8 -bit  Computer  1 nbs 
A-D/D-A  converter  sub-system,  a phase-locked  oscillator/synthesizer 
(PLO/SYN),  a Tektronix  4012  graphics  terminal,  a Versatec  1600A  high  speed- 
high  resolution  electrostatic  simultaneous  printer/plotter,  a Texas  Instruments 
Silent  700  terminal,  and  a Remex  high  speed  paper  tape  reader. 

The  original  design  goal  for  this  digital  signal  processing  system  was 
the'  spectrum  analysis  on  a continuous  real-time  basis  of  a 1.5  MHz  receiver 
fl'-HW  with  a Fast  Fourier  Transform  (FFT)  filter  frequency  resolution  of  300 
to  400  Hz.  The  first  design  therefore  included  a hardwired  array  processor 
with  a 1.5  million  word  throughput  rate.  Sponsor  induced  program  changes 
after  the  computer  portion  of  this  system  was  purchased  necessitated  systen, 
redesign  to  a multi-purpose  signal  processing  laboratory. 

The  present  system  will  digitize  analog  data  of  bandwidths  of  up  to 
100  KHz.  The  system  has  three  phases,  acquisition,  pulse  analysis  and 
plotting.  In  the  acquisition  phase  the  system  digitizes  analog  data,  packs 
the  data  samples  into  computer  words,  reads  the  data  into  computer 
memory,  unpacks  the  data  and  places  the  data  on  disk  memory.  During  the 

Note:  Manuscript  submitted  April  22,  1977. 
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pi-ise  analysis  phase  the  data  is  recalled  from  disk  as  needed  and  analyzed 
via  software  FFT  algorithms.  These  transforms  are  displayed  as  power 
versus  frequency  plots  on  the  graphics  terminal  and  also  stored  in  disk 
memory  if  desired.  In  the  final  phase  the  transforms  are  plotted. 

Hardware  Design 
A.  Analog-to-Digital  Controller 

The  major  in-house  hardware  design  effort  has  been  the  design  of  the 
Analog-to-Digital  controller  interface  circuits  (Fig.  3).  This  interface  unit 
contains  the  logic  for  controlling  the  A-D  converter  clock,  buffering  the 
3-bit  byte  data  from  the  A-D  converter,  generating  a byte  parity  bit,  and 
packing  the  8-bit  byte  data  sample  into  a 16-bit  computer  word  plus  byte 
parity.  The  interface  unit  then  sends  a request  to  send  to  the  computer. 

The  process  starts  with  an  end  of  conversion  or  Data  Ready  pulse  from 
the  A-D  converter . This  Data  Ready  pulse  is  reshaped  by  Schmitt  trigger  and 
one-shot  circuits  and  is  used  to  drive  a toggle  flip-flop.  The  outputs  of  this 
flip-flop  are  Anded  with  the  Data  Ready  pulse  to  provide  a shift-in  pulse  tc 
each  First  In-First  Out  Buffer  (FIFO)  in  turn,  starting  with  the  low  byte  FIFO. 
The  8-bit  data  sample  from  the  A-D  converter  is  sent  to  a parity  generator 
and  to  both  FIFO  buffer  inputs  . A parity  bit  is  added  to  the  data  sample  by 
the  parity  generator  and  the  selected  FIFO  buffer  is  pulsed  to  gate  the  data 
sample  and  parity  bits  in.  The  FIFO  buffers  are  9-bits  wide  and  32  words  deep. 
Two  of  these  units  are  paralleled  to  form  the  computer  word.  Data  gated  into 
either  FIFO  buffer  ripples  through  to  the  output  register  of  the  FIFO  and  sets 
the  output  ready  line  high.  This  output  ready  line  is  Anded  with  the  output 
ready  line  of  FIFO-2  so  that  as  soon  as  a word  is  formed  by  the  FIFO's  a re- 
quest to  send  is  sent  to  the  CPU.  After  acknowledging  the  request  the  CPU 
clocks  the  data  into  computer  memory. 

To  output  data  to  the  interface,  the  CPU  tests  the  state  of  the  output 
suffer  ready  line;  if  the  line  is  true  (meaning  the  buffer  is  empty)  the  data 
is  clocked  into  the  data  out  FIFO's.  The  input  ready  lines  of  these  FIFO's 
>»  are  Anded  and  the  output  of  the  And  gate  tied  to  the  above  sense  line.  The 

data  again  ripples  through  ooth  FIFO's  and  is  clocked  out  to  the  D-A  converts 
ov  an  external  clock.  Data  can  thus  be  readout  at  any  convenient  rate. 
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Fig.  3 — A-D  converter/computer  interface 
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FIFO  buffers  used  at  present  have  a 500K-word  minimum  thru-put  rate, 
so  that  the  maximum  A-D  conversion  rate  is  1 mega-samples  per  second. 

We  plan  to  replace  these  FIFO's  with  1 -mega -word-rate  units,  and  thereby 
increase  the  A-D  conversion  rate  to  2 -mega-samples  per  second.  The 
present  design  uses  a Non-Direct  Memory  Access  (DMA)  controller  and  is 
limited  to  the  Non-DMA  word  rate  of  225K-samples  per  second.  However, 
line  termination  problems  in  the  breadboard  proto -type  interface  now  in 
use  has  limited  the  present  usable  A-D  conversion  rate  to  about  half  this 
value . 

B.  Triggerable  Digital  Delay-Line 

The  Triggerable  Digital  Delay-Line  (TDDL)  is  a pulse  triggered  clock 
frequency  divider  with  a delay  range  from  1 milli-second  to  99.999  seconds, 
and  with  start-stop  outputs  to  drive  a time  interval  counter  (Fig.  4).  The 
stop  output  also  starts  the  A-D  converter  clock.  The  input  clock  frequency 
can  by  varied  from  1 MHZ  (tape  reference  frequency  at  12  0 IPS)  to  125  KHz 
(tape  reference  frequency  at  15  IPS).  Lower  clock  frequencies  (tape  speeds 
from  7\  IPS  to  1 3/4  IPS)  can  be  used  by  multiplying  the  TDDL  dial  settings 
by  the  factor  of  tape  speed  slowdown  less  than  15  IPS.  Delay  times  are 
synchronized  with  tape  speed  variations.  This  is  most  important  since  it 
allows  the  digitizer  to  be  started  at  the  same  event  on  tape  regardless  of 
tape  speed. 

TDDL  Circuit  Description 

A trigger  pulse  from  the  1 PPS  output  of  a time-code  reader  is  applied 
to  the  trigger  input  circuits  where  it  is  reshaped  and  sent  to  the  set  input 
of  a trigger  start-stop  latching  flip-flops.  The  Q output  of  this  latch  is 
Anded  with  the  clock  signal.  The  output  this  And  gate  is  sent  to  the  set 
input  of  the  clock  start  latch  whose  function  is  to  allow  only  the  first 
clock  pulse  after  triggering  to  be  sent  to  the  time  interval  start  one-shot. 
Next,  a binary  counter  receives  the  clock  pulse.  The  counter's  outputs 
are  each  sent  to  an  And  Gate.  The  other  input  of  each  Gate  is  tied  to 
a front  panel  switch  which  allows  one  to  select  the  proper  counter  output 
for  a 125  KHz  clock  output  signal.  This  clock  signal  is  divided  by  125  to 
give  a 1000  pps  signal.  The  1000  pps  signal  is  routed  to  the  next  divider 
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and  to  one  input  of  an  8-input  Nand  Gate.  Four  of  the  other  inputs  are 
from  later  stages  in  the  divider  chain.  The  sixth  input  is  the  output  of 
a SN7442  decimal  decoder.  These  five  inputs  are  normally  high  (at  a logic  one 
so  that  the  output  of  the  Nand  Gate  is  the  1000  pps  signal,  which  is  sent 
lO  a 7490  decade  counter.  The  decade  counter  outputs  are  decoded  and 
sent  to  a front  panel  switch.  This  switch  selects  the  appropriate  decoder 
output  for  the  desired  milli-second  count.  The  switch  output  is  sent 
to  the  above  mentioned  Nand  Gate,  inverted  and  sent  to  another  Nand 
Gate.  The  second  Nand  Gates  output  is  sent  to  the  time  interval  stop 
one-shot.  The  stop  one-shot  stops  the  time  interval  counter,  starts  the 
A-D  converter  clock, resets  the  counter/divider  chains,  and  resets  the 
trigger  start  latch  which  stops  the  clock.  The  other  circuits  are  identical 
and  give  selectable  counts  of  10  and  100  milli-seconds , and  counts  of  1 and 
10  second  steps.  The  entire  process  is  repeated  at  the  next  trigger  pulse. 
Normal  Operations  (acquisition) 

This  section  describes  procedures  to  scan,  locate  and  digitize 
taped  analog  signals.  The  analog  tapes  used  are  presumably  recorded  at 
120  ips  with  bandwidths  of  up  to  2 MHz.  Time  codes  and  a reference 
frequency  of  1 MHz  are  normally  recorded  at  the  same  time  as  the  date . The 
tape  track  of  interest  is  up-converted  to  10  MHz  and  applied  to  a bank  of 
100  analog  filters.  The  outputs  of  each  filter  is  detected  and  sampled,  with 
tne  sampled  output  displayed  on  an  oscilloscope.  The  times  and  types  of 
spectral  occurances  are  then  recorded.  The  exact  times  of  spectral  occurances 
are  found  by  recording  the  sampled  output  and  the  time  code  on  separate 
channels  of  a Bill-Howell  strip-recorder.  A 100  Hz  pulse  train  is  recorded 
on  another  channel  at  the  same  time.  Bandwidths  and  time  duration  of 
spectra  observed  is  noted  and  a decision  is  made  as  to  the  bandwidth  to 
be  used  for  digitizing.  Determining  factors  in  the  choice  of  the  bandwidth 
to  be  digitized  are:  spectra  bandwidth  and  time  duration,  initial  SNR 
translator  bandwidth  (40  KHz  at  present),  and  available  computer  memory. 

If  one  assumes  that  the  spectra  of  interest  is  centered  in  the  trans- 
lator passband  , the  output  of  the  translator  at  the  40KHz  6W  setting  is 
a signal  with  a bandwidth  of  from  4 to  44  KHz.  This  signal  is  low-passed 
at  44  Khz  and  applied  to  the  A-D  converter.  At  the  same  time,  the  time 
code  signal  from  the  tape  recorder  is  FM -Demodulated  and  sent  to  a time 
code  reader.  The  1 PPS  output  of  this  reader  is  used  to  trigger  the  TDDL. 


The  reference  signal  from  the  tape  is  filtered  and  sent  to  a phase-locked 
oscillator  synthesizer  (PLO-SYN)  and  to  the  TDDL  clock  input.  The  output 
of  the  PLO-SYN  is  the  A-D  converter  clock,  while  outputs  from  the  TDDL 
start  and  stop  a Hewlett-Packard  Time-Internal  counter  and  trigger  the 
gated  A-D  converter  clock  circuit. 

Once  the  appropriate  delay  is  decided  upon,  the  TDDL  is  set,  the 
acquisition  program  is  loaded  into  the  computer,  and  the  analog  tape 
recorder  is  started.  The  tape  recorder  is  normally  operated  at  a 8:1  speed 
reduction  (15  IPS)  and  one  has  time  to  arm  the  A-D  converter  interface  so 
that  tne  next  pulse  from  the  1PPS  time  code  reader  output  will  trigger  the 
TDDL  which,  in  turn,  will  start  the  A-D  converter  at  the  proper  time  interval. 
The  A-D  converter  interface  packs  the  8-bit  data  samples  into  16-bit  words, 
adds  the  byte  parity  bits  and  signals  the  computer  to  receive  the  data  samples. 
After  the  computer  memory  block  is  filled,  or  the  number  of  data  samples 
desired  taken,  the  computer  will  either  stop  the  A-D  converter  clock  or 
ignore  further  requests  to  send  data  from  the  A-D  converter  interface. 

The  tape  recorder  is  then  stopped  and  the  graphics  terminal  can  be 
observed  to  see  if  the  spectra  of  interest  was  digitized.  If  desired  the 
entire  process  can  be  repeated.  Normally  the  TDDL  settings  would  be 
adjusted  so  as  to  move  the  digitizer  time  window  to  an  adjacent  time  slot. 

Witn  repeated  passes  one  can  slide  the  digitizer  time  window  through  the 
time  duration  of  the  signal  of  interest. 

Alternate  means  of  data  acquisition  are  putting  the  signal  of  interest 
on  a video  disk  and  using  the  disk  sync  pulse  to  trigger  the  TDDL  or  if 
the  conversion  rate  is  low  enough,  continuous  digitizing  with  the  computer 
transferring  data  from  main  memory  to  disc  memory  and  to  magnetic  tape. 

Software  & Programming  Considerations 

The  Varian  V-73  mini-computer  system  (Fig.  5)  utilizes  a batch 
processing  operating  system  (N10S -version  K)  with  a Fortran  compiler, 
a DAS  assempler,  a text  editor,  and  a file  system.  We  have  expanded 
the  operating  system  to  include  the  Tektronix  Plot -10  software  package, 
and  a U.S.I.  nrdified  version  of  Versatec's  plotting  software.2 

Installation  of  vendor  supplied  software  and  the  subsequent  gen- 
eration of  routines  particularly  suited  to  our  signal  processing  application 
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Fig.  5 — System  components 


was  made  unduly  difficult  by  the  lack  of  a single  vendor  for  the  entire 
system.  Individual  software  packages  had  to  be  integrated  into  the 
Varian  operating  system.  In  particular  the  very  sophisticated  Versatec 
plotting  software  had  to  be  hand  tailored  for  the  Varian,  necessitating  a 
time  consuming  process  of  code  translation,  a mulitude  of  Fortran 
modifications  & careful  checking  of  test  routines.  Lengthy  negotiations 
with  Versatec  resulted  in  the  acquisition  of  a Varian  compatible  version 
of  the  Universal  Vesatec  Software  prepared  by  Underwater  Systems,  Inc. 

At  the  initial  stages  of  the  operation  it  was  necessary  to  generate 
assembly  language  code  to  select  an  input  gate  in  the  buffered  1/0  Con- 
troller and  sense  the  state  of  the  output  control  pulse.  A ready  state  re- 
sulted in  the  input  buffer  of  the  controller  being  read  into  memory.  The 
debugging  of  this  buffered  driver  occurred  simultaneously  with  the  check 
out  of  the  hardware  A-D  controller  interface  circuits.  Signal  transmission 
into  the  Varian  via  the  A-D,  and  back  out  via  the  D-A  was  monitored  with 
an  oscilloscope  to  test  the  accuracy  of  this  dual  hardware/software  effort. 

All  of  the  processing  was  completed  under  the  restraint  of  limited 
central  memory;  the  operating  system  generally  left  not  more  than  12  k of 
core  for  programs  and  data.  Therefore,  another  time  consuming  phase  of 
the  signal  processing  effort  was  the  efficient  use  of  the  disk  for  storage 
of  intermediate  results.  Numerous  utility  programs  were  written  yielding 
various  methods  of  storing  and  retrieving  the  data,  and  displaying  results 
alternately  on  the  graphics  terminal  or  the  line  printer. 

A recent  task  designed  to  demonstrate  the  capabilities  of  this  signal 
processing  laboratory  entailed  the  acceptance  of  data  from  an  analog  tape 
and  subsequent  analysis  and  display.  The  limited  amount  of  central  mem- 
ory available,  and  the  desire  to  process  long  strings  of  data,  necessitated 
a three  phase  operation,  with  the  data  being  stored  on  disk  in  the  interim. 

Phase  I (TEKTRAN)  consists  of  accepting  the  incoming  digitized 
signals,  unpacking  the  data  and  storing  the  raw  data  on  disk.  If  desired, 
the  data,  at  this  point,  can  be  biased  and  displayed  on  the  graphics  terminal. 
During  this  phase  the  Fourier  Transform  of  the  signal  is  computed  and  also 
stored  on  disk.  The  option  of  displaying  the  transform  on  the  graphics  term 
inal  is  also  available.  (Fig.  6) 
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Phase  II  (EIGTRAN)  consists  of  reading  sequential  files  from  disk 
(each  data  file  being  followed  by  its  transform)  and  creating  a plot  file 
for  each  data  file.  Storage  limitations  require  the  successive  readings 
of  small  amounts  of  data,  until  the  entire  data  set  is  scanned,  com- 
puting a scaling  factor,  and  then  repeating  the  entire  process  using  the 
computed  scaling  factor  to  generate  a plot  file  for  each  data  file  (Fig.  7). 

The  final  phase  involves  the  execution  of  VPLOT  (a  modified  Versatec 
routine  incorporated  into  the  MOS  operating  system)  to  convert  the  pen 
movements  delineated  in  the  plot  file  to  the  ordered  raster  output  which 
is  necessary  for  the  electrostatic  plotter.  The  plots  can  be  generated 
immediately  or  stored  for  output  at  a later  time  (Fig.  8). 


Conclusions 

A mini-computer  based  digital  signal  processing  laboratory  was 
discussed  in  terms  of  hardware  interface  design.  Pertinent  hardware 
and  sofware  design  and  performance  details  were  noted  and  a brief 
description  of  normal  operations  given.  While  the  present  signal  pro- 
cessing capability  is  limited  to  signals  of  100  KHz  EW  or  less,  the 
system  can  be  easily  expanded  by  the  addition  of  DMI  interface  cards 
to  process  signals  with  bandwidths  of  up  to  2 MHz. 


Varian  V-73  System  Handbook.  Varian  Data  Machines,  Irvine,  CA. 
Versaplot  Adaptation  to  VDM,  Underwater  System,  Inc. 
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Program  Listings 

PACE  1 

1 * TEKTRAN 

2 * 

3 * REVISION  I - VARYING  TITE  RESOLUTION  OF  THE  DATA  AND  THE 

4 * TRANSFORM 

»)  ^ 

6 * PROGRAM  ACCEPTS  DATA  FROM  A-D  CONVERTER  — 1024  PTS  PER  PASS 

7 * UNPACKS  THE  DATA 

0 * DISPLAYS  THE  DATA  ON  THE  TEKTRONIX 

9 * STORES  THE  DATA  ON  DISK  L.U.  14 

10  * COMPUTES  TRANSFORM 

11  * STORES  TRANSFORM  ON  DISK  L.U. 14 

1 2 * 

i:l  * 

14  * 

15  EXTERNAL  IN.DA.DISPLA 

I *.  1 NTF.GER  DATA2 

17  COMMON  DATA2( 1024) 

18  REWIND  14 

10  WRITE! 3. 88) 

20  88  FORMAT! 5411  THIS  PROGRAM  WILL  ACCEPT  DATA  FROM  THE  A-D  CONVERTER, 

21  1/.36H  COMPUTE  TRANSFORM,  AND  STORE  DATA  AND  TRANSFORM  ON  DISK/) 

22  * 

23  **********  EXTRACTS  N POINTS  ( POWT.R  OF  TWO)  FROM  DATA  SET********** 

24  * 

25  10  WRITE!  3.87) 

2»>  87  FORMAT! 33H  ENTER  NUMBER  OF  DATA  POINTS  TO  TRANSFORM  (13  FORMAT)/ 

27  X) 

28  READ! 3, 80 >N 

20  80  FORMAT'  13) 

30  IF(N.EO.O)  GO  TO  1000 

31  WRITE! 3.84) 

32  84  FORMAT!  66H  IF  LINE  PRINTER  OUTPUT  OF  DATA  IS  DESIRED.  TYPE  1 . IF 

33  X NOT.  TYPE  0) 

34  READ!  3.83)  LPFLAG 

35  85  FORMAT!  I 1 ) 

36  WRITE! 3. 83) 

37  83  FORMAT!  3411  IF  TEKTRONIX  OUTPUT  IS  DESIRED  TYPE  1.  IF  NOT,  TYPE  0 

38  X) 

39  READ!  3. 83)  ITF.KFL 

40  CALL  IN 

41  CALL  DA 

42  IF!  ITEKFL. EQ.O)  GO  TO  92 


i 


PACE 


o 


43 

42 

1 F AC- 1024/N 

44 

DO  100  JJ= 1.1024 

4.1 

JK=.IJ*IFAC-  IF  AO  l 

46 

100 

DATA2)  .JJ)  -DATA21  JK) 

47 

CALL  DlSPl.Al  N) 

48 

DO  2o  1 = 1 . N 

44 

20 

D ATA21  1 ) = DATA21  1 ) *8 

30 

IF(  LPFLAC.EQ.O)  GO  TO  101 

3 1 

MUTE!  3.40)  N 

32 

40 

1 ORMATf  IHI.4H  N=  . 15//) 

33 

WH 1 TE*  3,41)  ( DATA2!  K) . K= 1 . N) 

34 

4 1 

FORMAT)  1018) 

33 

10  l 

MUTE)  14)  ('  DATA21  K)  . K=  1 . N) 

36 

ENDriLE  14 

37 

C ALL  F'lUUK.  DATA2) 

38 

IF)  LPFLAC.EU.O)  GO  TO  102 

34 

WHITE) 5.40)  N 

60 

WR [ TE) 3 . 4 1 ) 1 DATA21 K) . K= 1 . N) 

6 1 

102 

WHITE)  14)  ) DATA2) K>  , K=  1 . N > 

62 

ENDFILE  14 

63 

CALL  SCAN) DATA2. DATA2. -N.350) 

64 

CALL  NODE)  -8.  XMIN.  DX.  XORG ) 

63 

CALL  N))DE(  -4.  YNIN.DY.  YORG) 

66 

WRITE) 3.78)  XMIN. YNIN. DY. YORG 

67 

78 

FORMAT* 2F30.4./.2F30.4) 

68 

1 DY- DY 

64 

IF) IDY.EO.O)  CO  TO  10 

70 

JMIN= YNIN 

71 

MAX*  lo-MDY+JMIN 

72 

1 1)Y=  ) MAX-JMIN)  /740 

73 

D))  30  I = 1 . N 

74 

30 

D \TA2<  I ) = < DATA2)  I)-JNIN)/IDY 

73 

IF)  LPFLAG.Ell.O)  CO  TO  103 

7 6 

WHITE) 5.40)  N 

77 

WH l TE) 3.41)  ( DATA2) K) . K= 1 . N) 

78 

103 

DO  31  I * 1 . N 

74 

31 

D \TA2)  1 ) = DATA2*  I) -300 

80 

IF)  ITT.RFL.  EO.  0)  GO  TO  10 

81 

CALL  DISPLA1N) 

82 

GO  TO  11) 

83 

! 000 

STOP 

84 

END 

1 lL» L. 


1 


1 * 
2 * 


SUBROUTINES  I N , DA,  I) ISPLA 


a * REVISION  I - VARIABLE  GRAPHIC  OUTPUT 

4 * 

5 * 

6 * 'IN’  \CCEPTS  DATA 

7 * 'DA'  UNPACKS  DATA 

» * ’D  ISPLA*  YIELDS  GRAPHIC  OiriTUT  ON  TEKTRONIX 

4 * 


10 

EXT 

TPL0.CH0U.s5SE 

1 1 

NAME 

DA 

12 

NAME 

IN 

13 

NAME 

D ISPLA 

000000 

c 

14  DATA2 

COMN 

1024 

000000 

000000 

A 

13  IN 

EN1R 

OO00O 1 

100060 

A 

16 

EXC 

0060 

000002 

ooooao 

A 

17 

LDXI 

-312 

0O0)>03 

1 77000 

A 

01(0004 

10 1060 

A 

ia  loop 

SEN 

0060. *+3 

1)00005 

OOOO 1 1 

R 

000006 

005000 

A 

10 

NOP 

OooooT 

00 1 000 

A 

2o 

JMP 

*-3 

OOOO  10 

000004 

R 

OOOO | | 

102060 

A 

21 

IMF. 

060, TEMP 

0000  1 2 

0o*I302 

R 

ooooi  a 

006027 

A 

22 

LDBE 

TEMP 

OOOO  14 

OO0302 

R 

000015 

006063 

A 

23 

STBE 

DATA+3 12, 1 

OOOO 1 6 

00 l 304 

R 

000017 

005 1 44 

A 

24 

IXR 

000020 

OO 1 046 

A 

23 

JXNZ 

LOOP 

00002 1 

000004 

R 

000022 

OO 100O 

A 

26 

RF.TU* 

IN 

00002a 

1 00000 

R 

000024 

000000 

A 

27  DA 

ENTR 

000025 

006030 

A 

2B 

LDXI 

-312 

000026 

177000 

A 

000027 

006020 

A 

20 

LDB I 

- 1024 

ooooao 

176000 

A 

0000a 1 

00740 1 

A 

30 

SOF 

000002 

0O6O  1 5 

A 

31  BACK 

LDAE 

DATA+3 12.1 

ooooaa 

O01304 

R 

000004 

00 1 00 1 

A 

32 

JOF 

RIGHT 

PACE 


1 


a 


000035 

00005 l 

R 

000o3o 

006  157 

A 

33 

LEFT 

ANAE 

MASKL 

000037 

oO 1304 

R 

000040 

004350 

A 

34 

LSRA 

0 1 0 

0«MI04  1 

006056 

A 

33 

STAE 

DATA2+ 1024 , 2 

000041! 

00200O 

C 

000043 

003 122 

A 

36 

IBR 

000044 

005 144 

A 

37 

IXR 

000045 

00  1046 

A 

311 

JXNZ 

DA+5 

000040 

00003 1 

R 

000047 

00 1 000 

A 

30 

JMP 

OUTPUT 

000050 

000060 

R 

00005 I 

006 157 

A 

40 

RIGHT 

ANAE 

MASKR 

000052 

001303 

R 

000053 

006056 

A 

41 

STAE 

DATA2+ 1 024 , 2 

000054 

002000 

C 

000055 

003 122 

A 

42 

IBR 

000050 

00 l OOO 

A 

43 

JMP 

BACK 

000057 

000032 

R 

OOOOoO 

006030 

A 

44 

OUTPUT 

LDXI 

~ 1024 

00006 1 

176000 

A 

45  * DATA  HAS  BEEN  UNPACKED  THIS  SECTION  MASKS  OUT  SICN 
BIT  ANDSHIFTS  IT. 


46 

* SETS 

UP  IX 

AND  IV  ADDS  BIAS  TO  IY 

THEN 

CALLS  TPLO 

000062 

0060 1 5 

A 

47 

LDAE 

DATA2+ 1024. 1 

000063 

002000 

C 

000o<)4 

Oo6  130 

A 

411 

ANA  I 

MASKS 

000065 

OO0200 

A 

O00O66 

oo425o 

A 

40 

LRLA 

0 MOVE  SIGN 

BIT 

INTO  PLACE 

000067 

OO 1010 

A 

30 

JAZ 

*+8 

000070 

0OO077 

R 

00007 l 

0060 15 

A 

5 1 

LDAE 

DATA2+ 1 024 . 1 

0ooo72 

002000 

C 

o0oo73 

006  1 10 

A 

32 

ORA  I 

0177400 

OOO-  >74 

177400 

A 

oooo75 

006055 

A 

53 

STAE 

DAT  12+  1 024 . 1 

000076 

0O2O00 

c 

OOOo77 

oo5 1 44 

A 

34 

IXR 

OOO 1 oO 

00 1046 

A 

55 

J XNZ 

OUTPUT-*- 2 

OOO |o 1 

000062 

R 

OOO 

OO 1000 

A 

36 

RF.TU* 

DA 

OOO  1 3 

1 00024 

R 

OOO | o4 

OoOOOO 

A 

57 

DISPLA 

FNTR 

Ooo l o5 

0o2o0o 

V 

511 

CALL 

USE 

21 


i 


i 

\ 

i 


H 


PACE 

a 

000106 

000000 

E 

000107 

OOoOO 1 

A 

50 

DATA 

1 

non  l m 

60  NP 

BSS 

1 

0OO  1 1 i 

006010 

A 

6 1 

LDA I 

1024 

<100  1 12 

0O2000 

A 

000  1 10 

147000 

1 

62 

SUB 

( NPl* 

<<00  1 14 

054166 

A 

63 

STA 

TEHP2 

000  1 15 

OO6030 

A 

64 

LDXI 

- 1024 

OOO  1 16 

176000 

V 

OOo 1 17 

006020 

A 

65  L00P2 

LDB1 

0 

000120 

OoOOOo 

A 

66  * CLEAR  SCREEN.  SET  ITYPE=0 

OOO 1 2 1 

006010 

A 

67 

LDA  I 

27 

000122 

ooooaa 

A 

000120 

054 1 56 

V 

68 

STA 

TEMP 

000124 

002000 

4 

6‘! 

CALL 

ClKlU.  TEMP 

OOO 1 25 

OO0O00 

E 

000126 

000002 

R 

000127 

OOoO 10 

A 

70 

LDA1 

12 

0001 00 

OOOo 14 

A 

OOO 1 0 1 

054  1 50 

A 

7 1 

STA 

TEMP 

000102 

002000 

A 

72 

CALL 

CHOU, TEMP 

OOO 100 

000125 

E 

OOO 1 04 

<100302 

R 

000105 

<1060  10 

A 

73 

LDA  I 

4 

OOO  1 06 

<100004 

A 

000107 

<154  142 

A 

74 

STA 

TF.flP 

000140 

<<<12<10<1 

A 

75 

CALL 

DELAY. TEMP 

000141 

00025 1 

R 

000142 

<100302 

R 

OOO 1 40 

00500 1 

A 

76 

T7.A 

000144 

006057 

A 

77 

STAE 

I TYPE 

OOO 1 45 

<10  1 306 

R 

000146 

00504 1 

A 

78 

TXA 

000147 

124133 

A 

70 

ADD 

TEMP2 

000150 

Oo 1 002 

A 

80 

JAP 

R 

OOO 1 5 1 

<100215 

R 

000152 

0060  15 

81  LOOP  1 

LDAE 

DAT \2+ 1 024 . 1 

000150 

00200O 

C 

000154 

006  I2i< 

A 

82 

ADD  I 

BIAS 

000155 

000606 

\ 

OOO  1 56 

<106057 

A 

83 

STAE 

IY 

OOO 1 57 

oo  1 3 1 2 

R 

r 


22 


PACK 


r 


\ 


out)  |60 

006067 

A 

84 

STBF. 

IX 

006 i 6 l 

(to  13  11 

K 

00O  1 02 

00502 1 

A 

83 

TBA 

ooo  i <>rt 

006 140 

A 

86 

SUB  1 

1023 

OOO  i<>4 

(to  1 777 

A 

(too  l<>5 

(to  1002 

A 

87 

JAP 

COPY 

000  166 

00o2  1 7 

R 

06616? 

O02000 

A 

88  PLOT 

CALL 

TPLO. I TYPE.  IX.  IY 

(too  17(1 

oooooo 

E 

00017 1 

00 1 306 

R 

000  1 72 

0O 1311 

R 

000  170 

001312 

R 

000174 

005144 

A 

80 

1XR 

OOO  175 

00504 1 

A 

00 

TXA 

000 170 

124104 

A 

0 1 

ADD 

TEMP2 

OOO 1 77 

00 | 002 

A 

02 

JAP 

COPY 

0(10200 

000217 

R 

00620 1 

00302 1 

A 

03 

TBA 

000202 

006 1 20 

A 

04 

AI)L»  I 

1 

000200 

00000 1 

A 

0002(14 

005012 

A 

05 

TAB 

OO0205 

00500 1 

A 

06 

T7  A 

000206 

003 1 1 1 

A 

07 

IAR 

000207 

0O6037 

A 

08 

STAE 

ITYPE 

000210 

00 1 306 

R 

0002  1 1 

OO l 046 

A 

00 

JXNZ 

LOOP  1 

0002  12 

Ooo 152 

R 

0002  1 0 

00 1 OOO 

A 

100 

JMP 

COPY 

0002 1 4 

0002  17 

R 

000213 

00 1 OOO 

A 

101  R 

RF.TU* 

DISPLA 

000216 

100104 

R 

0002  17 

0060 10 

A 

102  COPY 

LDA1 

105 

000220 

000131 

A 

00022 1 

034060 

A 

103 

STA 

TEMP 

000222 

002000 

A 

104 

CALL 

DELAY. TEMP 

000223 

0o023 1 

R 

000224 

000302 

R 

000225 

006010 

A 

103 

LDA I 

27 

000226 

(100033 

A 

000227 

054052 

A 

106 

STA 

TEMP 

000230 

002000 

A 

107 

CALL 

CHOU. TEMP 

00023 1 

OOO l 33 

F. 

000232 

(1(10302 

R 

p \gf. 


0002:13 

006010 

A 

toa 

LDAI 

23 

000234 

O(*0o27 

A 

000233 

034044 

A 

io<> 

STA 

TEMP 

000236 

002000 

A 

1 10 

CALL 

CHOU, TEMP 

000237 

O0O23 1 

E 

000240 

000302 

a 

00024 1 

0060 1 0 

A 

1 1 1 

LDA I 

105 

000242 

000  13  1 

A 

000243 

034036 

A 

1 12 

STA 

TEMP 

000244 

002000 

A 

1 13 

CALL 

DELAY, TEMP 

000243 

00023  1 

K 

000240 

000302 

a 

000247 

00  1 (lOo 

A 

1 14 

JMP 

LOOP2 

000230 

0001 17 

a 

00023  1 

00(1000 

A 

113  DELAY 

ENTR 

000232 

O0200O 

A 

1 16 

CALL 

SSE,  1 

000233 

ooO 106 

E 

000234 

OOOOO 1 

A 

000233 

1 17  TINE 

BSS 

1 

000236 

077000 

I 

1 1H 

STX 

SAVX 

000237 

067000 

l 

1 10 

STB 

SAVB 

000260 

006037 

A 

120 

LDXE 

( TIME)* 

00026 1 

100233 

R 

• >00262 

0060  1 o 

A 

121 

LDAI 

077777 

0002..3 

077777 

A 

1*00264 

0o33  1 1 

A 

122 

DAR 

(*00263 

003000 

A 

123 

NOP 

000266 

OO  10  10 

A 

124 

JAZ 

COUNT 

000267 

000272 

a 

000270 

00 1 000 

A 

125 

JMP 

*-4 

00027 1 

000264 

R 

000272 

003344 

A 

126  COUNT 

DXR 

(*00273 

00  1046 

A 

127 

JXNZ 

*-0 

000274 

000262 

R 

000273 

037000 

l 

128 

LDX 

SAVX 

000276 

027000 

I 

120 

LDB 

SAVB 

000277 

001000 

A 

130 

RETU* 

DELAY 

000300 

10023  1 

R 

00030 1 

00o777 

A 

131  X 

HIT 

0777 

000606 

A 

132  DIAS 

E(1U 

390 

000200 

A 

133  MASKS 

EOU 

12B 

000 1 77 

A 

134  MASKD 

EQU 

127 

000302 

133  TEMP 

BSS 

1 

000303 

136 

TEMP2 

BSS 

l 

000304 

137 

D vr  A 

ii>s 

3 12 

001304  177400  A 

1311 

NASKL 

D AT  A 

0 177400 

00  13o.">  00O377  A 

1 3) 

MASKR 

D \TA 

0377 

00 1306 

140 

1 TYPE 

BSS 

1 

001307 

14  1 

SAVX 

BSS 

1 

001310 

142 

SAvn 

RSS 

1 

001311 

143 

IX 

BSS 

l 

00 1312 
ENTRY  NAMES 

144 

143 

IY 

BSS 

END 

1 

000024  II  DA  000104  R DISPLA  000000  R IN 
EXTERNAL  NAMES 


♦100233 

SYMBOLS 

E 

sJSE 

000237 

E 

CHOU 

000253 

E 

SSL 

000032 

R 

BACK 

0002  1 7 

R 

COPY 

000272 

R 

COUNT 

000000 

C 

DATA2 

00025  1 

R 

DELAY 

001306 

R 

ITYPE 

00 1311 

R 

IX 

OOOO,.-; 

u 

LOOP 

000  1 52 

R 

LOOP  1 

00I3.J4 

R 

MASKL 

0!  1 ' 305 

R 

MAS  ICR 

00006 0 

R 

OUTPUT 

000 1 67 

R 

PLOT 

00  1 3 i 0 

R 

SAVR 

00 i 307 

R 

SAVX 

000253 

R 

TIME 

OOo 170 

E 

TPLO 

0 ERRORS  ASSEMBLY  COMPLETE 

0001 TO  E TPLO 


000606  A BIAS  000237  E CBOU 
000024  R DA  000304  R DATA 
000104  R DISPLA  000000  R IN 


001312  R IY 
000117  n LOOP2 
00O20O  \ MASKS 
00021;-.  R R 
000302  R TF.HP 
00o30i  R X 


000036  R LEFT 
000177  A MASKD 
1)00  110  R NP 
00005  1 R R l OUT 
000303  R TLMP2 
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t 
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4 


6 

■r 

» 

a 

9 

10 
1 1 
12 

13 

14 
13 
16 
17 
la 
i‘> 
20 
21 
22 

23 

24 
23 
20 
27 
2fl 
20 

30 

31 

32 

33 

34 
33 
30 

37 

38 

39 

40 
4 I 
42 


* BIGTRAN 

* 

***«:*. «:*:i:tf*********:ft**:K********W-S:***************S:*:;:**X*:t:***  ********* 

* HEAD  FROM  DISK 

* AND  PLOT  FROM  A FILE  - USING  SEGMENTS 

* FILES  ARE  OF  VARYING  LENGTHS.  KITH  THE  FIRST  FIEF.  CONTAINING 

* A FUNCTION  OR  DATA.  AND  HIE  FILE  IMMEDIATELY  FOLLOW  I NC  CONTAINING 

* ITS  FOURIER  TRANSFORM 

:R 

Hi 

■ii 

* VARIABLE  LENGTH  PLOTS 

* 

* REVISION  I:  FREQUENCY  SCALE 

* 

****tf:*:fc*:|c:|:*:fc*:f:**:l::t:**:t::l:**:fc:|c**:f::fc***tf*:K**tt:t:*:f:*:f:*:fc***rc****:|t****!K***:K:fc*****:(:* 


INTEGER  BUF 

DIMENSION  BUF(faO) . X! 00) . I ALF 1(8) . IALF317) . IALF413) 

D AT V SAME/9999./. I ALF2/2HY  /.  IALF0/2HX  / 

D \T A ( IALF1C  I)  . 1=  1 .a>/2HFR.2!IEQ.2HUF..2HN0.2HY  , 2IK  K.2HII7.2H)  / 
DATA  ( I ALF3I  I)  , 1=  1 . 7 > /2H0R.  211 1C.  2H  I N . 2HAL . 211  D.21IAT.2HA  / 

DATA  ( I ALF4(  I).  I = 1 . 3) /2HTR.  2IIAN . 2IISF . 2II0R.  2HM  / 

I PLOT= 1 
REWIND  14 
WRITE! 3. 80) 

30  FORMAT!  3.3H  PROGRAM  PLOTS  FUNCTION  AND  TRANSFORM  STORED  ON  L.U.  1 

X4/> 

300  WRITE.  3.86) 

06  FORMAT!  3 1 II  TYPE  PTS.  IN  FILE.  SAMPLING  FREQUENCY  (13.F10.0)  ) 

READ. 3 . 89 ) N . SAMFR 
89  FORMAT!  1 3. FI 0.0) 

* IF  N = O , NO  MORE  PLOTS  ARE  DESIRED 
IFiN.EU.O)  GO  TO  700 

WRITE! 3.84) 

84  FORMAT!  0.311  IF  LINE  PRINTER  OUTPUT  OF  DATA  IS  DESIRED  TYPE  1,  IF 

KNOT.  TYPE  0) 

READ* 3. 83 > LPFLAG 
83  FORMAT!  I I > 

* WHEN  THE  DATA  AND  TRANSFORM  ARE  STORED  ON  DISK  THE  SYSTEM 

* CRE.iTES  60  WORD  RECORDS  - THE  FOI.LOWTNC  CODING  DETERMINES 
HU.  NUMBER  OK  RECORDS  ( NREC)  CREATED 


26 


P \GF. 


o 


43  NUM=  N/60 

44  Mfl-  N-NUM*60 

45  M=60-MM 

46  I F( H)  100.101,100 

47  100  NUM= NUM+ 1 

40  lOl  NREC=NUM 

4‘>  L£N=N/50+l 

50  FLEN=  LEN 

51  I F!  N . LE. 350)  FLEN=7. 

52  CALL  MODE!  2. FLEN. SAME. SANE) 

53  * FREQUENCY  RANGE  FOR  TRANSFORM  EQUALS  ONE  HALF  THE  SAMPLING 

54  * FREQUENCY  OR  'FOLDING  FREQUENCY’ 

55  F0I.DF  = . 3*SAMFR 

56  FINPPT=FLEN/FLOAT! N-l) 

57  * JK= 1 IMPLIES  PLOT  OF  FUNCTION 

50  * JK=2  I MPL I ES  PLOT  OF  TR\NSFORN 

50  DO  600  JK= 1.2 

60  IF(JK. Eft.  1)  DELTAX- FLOAT!  N- 1 ) /FLEN 

61  IF! JK. EQ. 2)  DELTAX= FOLDF/FLEN 

62  * J=l.  SCANNING  OF  DATA  FOR  SCALING 

63  * J= 2.  PLOTTING  OF  SEGMENTS 

64  DO  500  J= I . 2 

63  IFLAG=J-1 

66  I F ( LPFLAG. Eft. 0)  GOTO  170 

67  1.11  ITF.t  3.00)  N.J.JK 

60  00  FORMAT!  1 II 1 . 4HN  = .13,311  .1  = .I2.6H  JK  = .12,/) 

60  170  NFS= I PLOT- 1 

70  I F( NFS . EQ. 0)  GO  TO  101 

71  DO  100  JF-l.NFS 

72  100  CALL  SKIP 

73  * LOCATES  DATA  FOR  NEXT  SEQUENTIAL  PLOT 

74  18 1 F K l = 0 . 

73  KN=60 

76  DO  200  1=1. NREC 

77  IF! I. EQ. NREC)  KN=-MM 

70  READ! 14)! BUF! K) , K= 1 . 60) 

70  IF! LPFLAC. EQ. 0)  GOTO  104 

00  WRITE!  3.9  1)  ! BUF! K) . K= 1 .60) 

01  01  FORMAT! 1010) 

82  104  DO  100  11=1.60 

83  X! I I)=FK1*FINPPT 

84  FK1=FK1+1. 


u;f. 
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ENTRY COMMON 


CONT I N1JE 

IF*  IF LAO)  198.197,198 
C U.I.  MODE!  7 . FLF.N . SAW..  SAME) 

C U.L  SCAN! X. BUF , KN . 050 ) 

THIS  CALL  SEQUENTIALLY  SCANS  BLOCKS  OF  DATA  — 

AFTER  LOOKING  AT  LAST  BLOCK.  SCALINC  FACTORS  ARE  DETERMINED. 
CO  TO  200 

0 \LL  MODE! 8.0. . DEL TAX. SAME) 

IF<  I.EQ.NRF.C)  GO  TO  196 
I F < I.EQ.  1 ) GO  TO  195 
CALL  DRAW!  X, BUF, 60. 039) 

GO  TO  200 

IN  ORDER  TO  ELIMINATE  A GAP  BETVEEN  PLOT  SECMENTS 
SUBSEQUENT  CALLS  TO  DRAW  ARE  WITH  THE  PEN  DOWN 
C VI. L I)RVW(  X.  BUF . 60 .051) 

THIS  CALL  SEQUENTIALLY  PLOTS  BLOCKS  OF  DATA 
Go  TO  200 

0 U.L  DRAWi  X. BUF. MM. 059, 

CONTINUE 

REWIND  14 
CONTINUE 

IF(JK.EQ.l)  CALL  AXES* 2.1.1 ALFO .2.1.1 ALF2) 

IF,  JK.F.Q.2)  CALL  AXES,  1 6 . 3 . I ALF  1 . 2 . 1 . I ALF2) 

CALL  MODE,  6. 4. SAME. SAME, 

CALL  MODE,  4 . . 2 , . 15. SAME) 

IFlJK.EQ.2)  GO  TO  502 
C \LL  NOTF.(  1.75,8.0,1  ALF3 .14) 

GO  TO  503 

CALL  NOTE,  2.  1.8.0,  I ALF4.  10) 

CALL  DRAWN  0.  , 0.  . I . 9000) 

THIS  CALL  ENDS  A PARTICULAR  PLOT 
WRITE, 3.87,  [PLOT 

FORMAT,  10K  PLOT  NO.  . I5.9H  COMPLETE) 

1 PLOT1  I PI.OT+  1 
RF.WIND  14 

GO  TO  300 
CONTINUE 

CALL  DRAW,  0.0. 0,999) 

THIS  CALL  INDICATES  THAT  ALL  PLOTTINC  IS  COMPLETED! 

STOP 

END 

BLOCK  NAMES 


» 


